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speckle  Interferometry 
For  a  Partially  Coherent  Source 


1.  INTRODUCTION 

The  conventional  theoretical  result^"®  for  the  spatial  frequency  spectrum  of 
a  source,  as  obtained  using  speckle  interferometry^’  ^  is  based  on  the  assumption 
that  the  source  is  spatially  incoherent.  This  assumption  is  not  restrictive  when 
imaging  distant  stars,  because  the  coherence  patches  are  generally  unresolvable 
by  the  imaging  system:  however,  when  one  tries  to  obtain  detail  over  sizes  on  the 
order  of  the  spatial  coherence  length,  the  conventional  theory  is  no  longer  valid, 
and  must  be  modified.  In  this  paper  we  will  obtain  the  frequency  transfer  function 
for  the  combination  of  atmosphere  and  imaging  system  in  the  case  when  the  source 
Is  partially  coherent,  and  can  be  represented  by  the  quasi -stationary  model  em¬ 
ployed  by  Carter  and  Wolf,  and  Leader.  “  The  limiting  case  when  the  source 

consists  of  2  point  sources  with  random  phase  relationship  has  been  studied  by 
12 

Miller  and  Korff. 


(Received  for  publication  29  March  1979) 

Due  to  the  number  of  references  to  be  Included  as  footnotes  on  this  page,  the 
reader  is  referred  to  the  list  of  references,  page  15. 
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2.  ANALYTICAL  FORMALISM 


6  7 

In  the  Labeyrie  *  procedure  a  series  of  short  exposure  photographs  of  a 
star  are  taken  and  then  successively  projected,  in  the  conventional  manner,  by  a 
laser  onto  a  second  piece  of  film.  This  leads  to  an  evaluation  of  the  ensemble 
average  of  the  modulus  squared  of  the  Fourier  transform  of  the  irradiance,  and 
gives  finer  detail  than  conventional  seeing  calculations  predict  should  be  resolv¬ 
able  in  the  presence  of  atmospheric  turbulence. 

In  order  to  discuss  the  aforementioned  method  quantitatively,  let  us  refer  to 

Figure  1.  If  the  field  distribution  of  the  source  is  denoted  by  u  (a,),  the  field  in 

o-*-l 

the  image  plane,  where  film  is  located  is, 
g-ik(L+x)  i?/*  „ 

“f<£>  JJ  (£v£^  (1) 

-flO 

where 

00 

^  d^Pg  T<£2)  exp  j  ik£2-  +'''<e2*£l>|  •  <2) 

In  writing  Eq.  (1)  we  have  assumed  that*  the  film  is  in  the  Fraunhofer  zone  of  the 
source,  defined  as  the  transmissivity  of  the  aperture  and  k  =  2»/X  as  the 
signal  wavenumber.  Also  ^^'e  additional  complex  phase,  due  to  atmos¬ 

pheric  turbulence,  of  a  spherical  wave  propagating  ft  oin  the  point  (x,£^)  on  the 
source  to  the  point  (0.jB2^  the  aperture.  The  short -exposure**  transparency  on 
the  film  is  proportional  to  N  such  short  exposures  are  made,  with  the  time 

interval  between  successive  exposures  greater  than  the  atmospheric  coherence 
lime,  so  that  we  have  N  statistically  independent  short  exposure  photographs.  If 
a  photographic  plate  is  then  placed  at  x  =  -L,  and  these  exposures  successively 
projected  onto  it.  It  is  found  that  the  transparency  on  the  plate  is  proportional  to 
<  |ij(l^/L)|^),  where 


If  the  film  is  in  the  Fresnel  zone,  we  simply  replace  u-(p)  by  u,(p)  exp  (lkp^/2L). 

I  ft 

By  short  exposure  we  mean  short  compared  with  the  coherence  time  of  the 
atmosphere. 
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I  LENS  WITH 


PLANE  A 


photographic  PLATE 


Figure  1.  Geometry  Assumed  for  the  Imaging  System 


JJ  I  Uj(£)  I  ^  e'ii  £■ 


(3) 


u  is  the  radian  spatial  frequency,  and  (  )  denotes  an  ensemble  average  over  the 
statistics  of  the  turbulence.  This  average  occurs  because  the  transparency  on 
the  plate  is  the  sum  of  N  »  1  short  exposures. 

We  can  relate  the  quantity  in  Eq.  (3)  to  the  source  distribution  by  substituting 
Eq.  (1)  for  Uj(£).  The  result  is 


Ijiu)  =  (2»Xx) 


-2 


u  L 


/•••/  d^pj  d2p2  d2p3  u^^^lu^^g)  T(£2)T-(£2 


exp  I  0(£2.^i)  +  <^'*(£2  ‘iL/k,£3)  +  i  -^£2  •  (£j  -£3)  -  i  •  £3^ 


(4) 


The  overbar  on  the  quantity  '“(,^£1^ ^£3^  indicates  an  average  over  source  statis¬ 
tics,  and  arises  because  it  has  been  tacitly  assumed  that  the  exposure  time  of  the 
individual  speckle  pattern  is  much  longer  than  the  source  coherence  time,  but  of 
course,  much  shorter  than  the  coherence  time  of  the  turbulent  atmosphere. 

For  a  spatially  incoherent  source  one  usually  assumes  ^q(£j)^q(£3)  = 

^£j  "£3^  where  6(. . .)  is  the  Dirac  delta  function  and  1^  is  the  source 
irradiance  distribution.  However,  for  partially  coherent  sources  this  assumption 
is  invalid,  and  instead  we  shall  employ  the  quasi-stationary  model  used  by  Carter 
and  Wolf,  and  others.  That  is,  we  shall  write 


V£i><i£3'  *  Is  8<£l  -£3> 

where  g  is  the  correlation  function  of  the  source.  In  writing  Eq.  (5),  it  is  tacitly 
assumed  that  the  source  coherence  length  is  much  smaller  than  the  total  size  of 
the  source. 
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We  now  use  Eq.  (5)  In  Eq.  (4),  define  sum  and  difference  coordinates 
^  *£l  '£3*  *  ^£1  ^£3V'2.  «<«■•,  W'd  perform*  the  ensemble  average  over 

the  turbulence  statistics,  assuming  that  the  source  lies  wholly  within  an  isoplanatic 
turbulence  patch,  so  that  ~  i4<p,0).  The  result  is 


H(^) 


(6) 


where 


OD 

i^(n)  •  (2»)*^  JJ  d^v  I^(v)  exp  (-in  •  v) 


(7) 


Therefore,  the  modulus  squared  of  the  Fourier  transform  of  the  source  distribu¬ 
tion  is  linearly  related  to  the  measured  quantity,  Also  in  Eq.  (t>),  we  have 

defined 


H(u))  =  (Xx) 


Jj  d^y  exp  | 

(i-HT-o)  *i  Di  (i'hf-  0)1  • 


where 


(8) 


90 

Ko(w,  1)  ’  //  d^T  T<T  ♦  2./2)  T*  (t  +  2./2  +  W  L/k) 

-oe 

•  T*<t  -  X/2)T(T  -  x/2  +  £L/k)G(T,2,,i£)G*(^  -2,u;)  ,  (?) 


*The  derivation  of  the  average, 

<exp(t(/^2*£i'  ^  '^'*<£5»£4'  ^  *<(£5  +iJL/k,£^)l> 

over  the  statistics  of  the  turbulence  is  given  by  Eq.  (34)  of  R.  Fante,  "Some 
results  on  the  imaging  of  incoherent  sources  through  turbulence,"  J.  Opt.  Soc. 
Am.  ^  574-580  (1976). 
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and 


Giry.u)  » 


oo 


?  Rif)  exp 


(T 


■>/2  ^  uL  2k)  I 


(10) 


1  lie  structure  function  Dj  of  the  atinosptieric  turbulence  is  Riven  bv 

D,(fi0)^  2|3/r^p/3  , 


(11) 


and  the  atmospheric  coherence  length,  r^,  is  defined  as 


r  = 
o 


(12) 


2 

where  Cj^(x')  is  the  index -of-refract ion  structure  constant  of  the  atmosphere. 

The  result  m  Eq.  (6)  is  of  the  same  form  as  previous  results  obtained  for  a 
spatially  incoherent  source,  but  differs  in  the  definition  of  K  (i£.x).  For  a  spa¬ 
tially  incoherent  source,  so  that  g(£j  *  £3)  =  6  -  £^)  where  is  the  co¬ 

herence  area  of  the  source,  we  find  G  «  A^,.  and  Eqs.  (6)-<10)  then  reduce  iden¬ 
tically  to  the  result  of  Korff  for  the  spatially  incoherent  source. 


3.  APPROXIMATE  EVALUATION  OF  TRANSFER  FUNCTION,  H(iC  ) 

We  would  now  like  to  investigate  the  effect  of  the  partial  coherence  on  the 
transfer  hinction,  H(u).  In  order  to  do  this,  we  will  assume  that  the  source  co 
herence  function  g(5)  is  given  by 


where  a  is  the  source  coherence  length.  Upon  using  this  expression  in  Eq.  (10), 
we  find 


>,  w> 


Ap  exp 


where  A_ 
c 


*a  is  the  coherence  area  of  the  source. 


(14) 
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The  next  step  is  to  substitute  G  into  Eq.  (9)  and  evaluate  K^(u,  7).  For  a 
realistic  aperture,  such  as  on  a  telescope,  "!  ‘.£2^  "  ^  ^£-2^  “  ^ 

for  Ipgl  >  D.  Unfortunately,  the  expression  for  K  in  Eq.  (9)  cannot  be  evaluated 
in  closed  form  when  this  expression  is  used  for  T(£2)-  Therefore,  we  have  been 
forced  to  approximate  T(£2)  by  the  function  reali¬ 

zation  that,  because  we  are  approximating  the  edge  diffraction,  our  results  will 
be  valid  only  for  radian  spatial  frequencies  such  that  |ij|  <  kD/L.  That  is,  setting 
T(£  )  =  0  for  |£2i  >  D  would  give  H(u)  =  0  for  |w|  >  kD/L:  however,  because  we 
will  approximate  T(£2)  by  exp  (-  we  will  obtain  an  expression  for  H(u) 

which  does  not  vanish  for  lw|  >  kD/L.  If  the  aforementioned  approximation  is 
used,  along  with  Eq.  (14),  in  Eq.  (9)  we  find  that 


where  f  =  2^/^x/(ka)  =  (2/»A^)^'^^  Xx  is  the  coherence  length,  as  measured  at 
the  receiving  aperture,  of  the  field  radiated  by  a  coherence  patch  of  radius  a  on 
the  source. 

We  next  investigate  the  behavior  of  the  Ainction 

C(w,7)  =  exp{-  Dj(uL/k,0)  -  0^(7,  0)  + -  uL/k,  0)  • 


It  is  found  that  for  L'l  «  kr  /L,  this  function  can  be  approximated  by  unity.  How 
ever,  the  more  interesting  case  occurs  when  »  kr^/L,  In  this  limit, 

C(u,  7)  can  be  approximated  by 

C(u,x)^«* **P  {■  2|x/*'ol^^^^  • 

The  result  in  Eq.  (17)  is  inconvenient  for  analytical  purposes;  consequently,  we 
shall  approximate  the  exponent  in  Eq.  (17)  by  a  quadratic  function,  so  that 


*The  evaluation  can  be  performed  In  the  limit  when  G(C)  =  A^. 

**We  are  assuming  that  ro  <  D.  For  r  »  D,  the  atmospheric  turbulence  has  an 
insignificant  effect  on  the  Imaging. 
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■)  '  ^  f  xp 


I  -  2.  :*'•(!  h  '\,1  ■  I  .  (i8> 

The  luni'tiun  in  Kq.  (18)  is  siu-h  that  equals  e  ^  at  the  same  value  of  111  as  does 
the  function  in  Kq.  (17). 

If  we  ni'w  use  Kqs.  (lii)  and  (18)  in  Kq.  (8),  we  can  calculate  the  transfer 
function  H(^),  {'or  values  of  |j£  1  such  that  |^|  S.S.  kr^^  1.,  the  result  is* 


(ly) 


(-’0) 


Hecause  of  our  approximation  for  T^.,).  the  result  in  Kq.  (IP)  can  be  applied  to 
realistic  apertures  only  for  spatial  frequencies  such  that  LI  v  kO/L.  In  the  limit 
when  r^^  v«.  n  and  the  source  is  incoherent,  so  that**  -  jc.  Kq.  (19)  reduces  to 
the  high-spatial  frequency  approximation  obtained  by  Korff.  ^ 

Kor  values  of  LI  such  that  y  v.  kr^^  1  the  evaluation  of  Kq.  (8)  gives^ 


(21) 


A  qualitative  plot  of  MU')  is  shown  in  Figure  2.  for  the  limits  when  r  w  O 
and  r^^  s.s  i>_  In  this  latter  case,  of  course,  the  turbulent  atmosphere  does  not 
affect  the  imaging. 


Note  that  HL')  is  not  normalized  to  unitv  at  u'  =0,  as  is  often  done. 

• 

Strictly  ?ipeaking  never  really  becomes  infinite,  because  a  '*  X.  Therefore 
f  2^  2  ,-l  X. 

'Because  r^,  is  assumed  to  be  leas  than  D,  the  assumption  that  LI  «*-  kr  /I. 
implies  that  L  L  '‘^^l  C'cnsequentlv,  in  B.'q.  (21)  we  have  approximated 

exp(-u;-l.2  by  unity. 
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gure  2.  Qualitative  Plot  of  H(u).  The  point  P  corresponds  to 
~  kr  /L  ~ 

—  o' 


4.  DISCUSSION 

From  Figure  2.  we  observe  that  the  amplitude  of  H(u)  at  small  spatial  fre¬ 
quencies  is  determined  by  the  size  of  D/f  For  f  ^  »  D  we  get  the  usual 
Labeyrie  result,  but  when  f  «  D  the  result  in  Eq.  (21)  becomes  H(u)  -  1. 

It  is  also  interesting  to  observe  the  nature  of  the  high -spatial-frequency 
plateau  in  H(ij)  for  the  case  when  r^  «  D.  We  will  not  discuss  the  case  when 
»  D  because  in  that  limit  atmospheric  turbulence  does  not  affect  the  image. 

If  r  «  D  but  f  »  D,  we  get  the  previous  result  of  Korff;  however,  if  r  «  D  and 
II  o 

t  «  D  we  find,  for  |u  >  kr  /L 


H(u,)  ^  (-u)=^L^/kV) 


1  +  2.296 


(22) 
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